Superconducting Single-Electron Transistor in a Locally Tunable Electromagnetic 
Environment: Dissipation and Charge Fluctuations 
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We have developed a novel system consisting of a superconducting single-electron transistor (S- 
SET) coupled to a two-dimensional electron gas (2DEG), for which the dissipation can be tuned 
in the immediate vicinity of the S-SET. Within linear response, the S-SET conductance varies 
nonmonotonically with increasing 2DEG impedance. We find good agreement between our experi- 
mental results and a model incorporating electromagnetic fluctuations in both the S-SET leads and 
the 2DEG, as well as low-frequency switching of the S-SET offset charge. 
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Electrical transport in nanoscale devices is strongly af- 
fected by the electromagnetic properties of their environ- 
ment. This is particularly true for superconducting sys- 
tems, and has been a topic of considerable interest lately: 
dissipation can drive a superconductor-insulator quan- 
tum phase transition jl], |^, |j| , and is also expected to af- 
fect the coherence time of qubits such as the charge-based 
single Cooper pair box [Q, |5|, 0]. Maintaining quantum 
coherence in such devices long enough to allow many op- 
erations is prerequisite for their use in quantum computa- 
tion. In this regard, the closely related superconducting 
single electron transistor (S-SET) is an excellent system 
for attaining a better understanding of the effects of the 
electromagnetic environment on quantum coherence. 

This approach was followed by the Berkeley group [Q , 
who fabricated an S-SET above a two-dimensional elec- 
tron gas (2DEG) in a GaAs/ALjGai-xAs heterostruc- 
ture. Using a back gate to vary the 2DEG sheet resis- 
tance R sq they measured changes in the S-SET conduc- 
tance Gset as the dissipation was varied. Building on 
earlier theoretical work || [| [ll], [l^], Wilhelm, et al. 
predicted []l3f that within linear response Gset would 
scale with the ground plane conductance G2D = 1/Rsq 
and temperature T as G^ D /T Q and that the S-SET cur- 
rent / at fixed bias voltage V would vary nonmonoton- 
ically with G2D, while switching from the non-linear to 
linear regime. These predictions were made within two 
models: in one the environment was treated as a ground 
plane, while in the other it was treated as an infinite 
RC transmission line provided by the SET leads. While 
the Berkeley group did observe power law behavior, their 
measured exponents were not in quantitative agreement 
with theory. Furthermore, the measured (3 depended on 
T and a on G2D , calling the scaling form into question. 

In this Letter, we report measurements on samples sim- 
ilar to those studied by the Berkeley group; in our sam- 
ples, however, we can modify G2D in the immediate vicin- 
ity of the S-SET while leaving the 2DEG beneath the 
leads largely untouched. In contrast to the predictions of 
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FIG. 1: (a) Electron micrograph showing the S-SET, gold 
gates, and QPCs. (b) S-SET I-V characteristics for V g = 0. 
The linear regime extends to roughly ±8 fiV. (c) Circuit 
diagram of the S-SET, including its environment. We allow 
both for lead impedances Zg and a ground impedance Z2D 
coupled to the S-SET through a capacitance C^d- 



Wilhelm, et al. we observe a nonmonotonic dependence 
of Gset on G2D entirely within the linear regime. We 
propose a model for the environmental impedance of S- 
SET/2DEG systems that includes electromagnetic fluc- 
tuations in both the 2DEG and leads, while treating the 
latter as finite RC transmission lines. We also allow for 
relatively low-frequency switching of the charge state of 
the SET island ^M, which can affect the measured cur- 
rent. Within this model, we find good agreement between 
our calculated and measured results. 

An electron micrograph of a typical sample is shown 
in Fig. |l](a). We begin with an GaAs/ALjGai-^As het- 
erostructure grown on a GaAs substrate using molecular 
beam epitaxy, consisting of the following layers: 1000 nm 
of GaAs, 47 nm of Alo.3Gao.7As and 5 nm of GaAs. 
The Alo.3Gao.7As is delta-doped with Si 22 nm from 
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the lower GaAs/Alo.3Gao.7As interface, at which forms 
a two-dimensional electron gas (2DEG) with R sq = 20 Q 
and sheet density n s — 3.6 x 10 11 cm -2 . On the sam- 
ple surface we use electron-beam lithography and shadow 
evaporation to fabricate an A1/A10 S-SET surrounded 
by six Au gates When no gate voltage is applied 

and the 2DEG is unconfincd, the measured I-V charac- 
teristics are linear over several microvolts, as shown in 
Fig. 0(b). We can also apply a single gate voltage V g to 
any combination of Au gates, excluding the 2DEG be- 
neath them. We focus on two geometries: the "pool," 
in which all six gates are energized, and the "stripe" in 
which only the four exterior gates are. In both cases, 
electrons immediately beneath the SET are coupled to 
ground by quantum point contacts (QPCs) with conduc- 
tances l/i?Qpc (assumed equal) as low as 3 conductance 
quanta Go = e 2 /h. In the stripe geometry, the electrons 
can also move vertically through a resistance i? s t r to a 
large 2DEG reservoir that is coupled to ground through 
a capacitance G str . As illustrated in Fig. 0(c), electro- 
magnetic fluctuations in the environment can couple to 
the tunneling electrons in two ways: through the leads, 
which act as RC transmission lines with impedance Zg for 
the relevant frequency range J?], [k| , and through the ca- 
pacitance C2D to the 2DEG with impedance Z^d, which 
is related to -Rqpc and (for the stripe) i? str . The model 
has been studied previously |l6|, [ijj without considering 
particular forms for and Zi. 

Measurements were performed on two separate sam- 
ples (SI and S2) in a dilution refrigerator in a four-probe 
voltage biased configuration; the estimated electron tem- 
perature was 100 mK. High frequency noise was ex- 
cluded using standard techniques. A small capacitance 
C g « 20.3 aF (not shown in Fig. 0(c)) couples the six Au 
gates to the S-SET. The other sample parameters such as 
the junction resistances and capacitances Gi^, the 
coupling capacitance G2D and superconducting gap A are 
given elsewhere [[t8|. The charging energy E c = e 2 /2Cs 
for sample SI (S2) is 118 (77) /xeV while the Josephson 
energy Ej. — A averaged for the two junctions is 

4.7 (21.8) /LteV. Here Rq = is the superconducting 
resistance quantum and Gs = C\ + G2 + G2D + C g . We 
use standard lock-in techniques and voltage biases of 3 
and 5 jj,V rms respectively to measure Gset and the con- 
ductance G2D across the series combination of the QPCs 
versus V g in the pool and stripe geometries. The results 
for S2 are shown in Fig. 

From Fig. ||(b), we see that for the pool Gset rises by 
nearly a factor of 2 as V g becomes more negative, before 
dropping rapidly. Although G2D vs. V g is nearly identi- 
cal in both cases, Gset for the stripe rises only by ~50% 
and does not decrease, even for the most negative V g . In 
both plots, e-periodic Coulomb blockade oscillations are 
seen as the S-SET offset charge varies. We fit a smoothly 
varying function to the measured G2D versus V g (inset, 
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FIG. 2: (a) Maxima of Gset vs. G2D in the stripe (o) and 
pool (A) geometries. Solid line: Gg ET vs. G2D- Inset: G2D 
(solid line) for the pool, and a c h in the stripe (x) and pool 
(O) geometries vs. V g . (b) Gset vs. V g in the pool geometry, 
(c) Gset vs. V g in the stripe geometry. 



Fig. g(a)) which we use to plot the maxima of Gset ver- 
sus G 2D in Fig. |(a). The knee in G 2D at V g « -0.38 V 
corresponds to the appearance of quantized plateaux in 
the individual QPC conductances. 

To understand these results and the I-V curve in Fig. 
(b), we begin with the rate of sequential Cooper pair 
tunneling fll9|] through junction j, valid for Ej <C E c 



T(5fW) = ^/2n)E%P(-Sf^) 



(1) 



where 6f® = // - /; = (~l) j 4E c (2N -n g + l)- 2a je V 
is the free energy difference for changing the num- 
ber of Cooper pairs N by 1, ctj = 5 + (— l) 3 (Gi — 
G 2 )/2G S , and n g = C g V g /e. Here P(E) is the prob- 
ability of exchanging energy E with the environment 
and can be expressed in terms of a correlation function 
K{t) = Rq 1 JZo ^Re[^ t ^)]{coth( 5 |^)[cos(^) - 1] - 
i sin(t^)} via P(E) = ^ dt exp[K(t) + i^] where 
Z t (uj) is the total impedance seen by tunneling electrons. 

The general result for Z t (uj) within the model of 
Fig. 0(c) is quite complex In our case, however, 

Zi (Z21)) dominates at low (high) frequencies and to an 
excellent approximation 



Be[Z t (w)] = /eiRe[Z*(w)] 



where for junction 1(2) ki 
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line Zrc(u>) — yj j^r t ta.nh(\/iuirgce£ 2 ) where m and eg 

are the resistance and capacitance per unit length and I 
is the total line length. We are interested in the long-time 
limit of K{t) which is dominated by the low- frequency 
part of Zt(u>). In that limit, Zrc{u) 
which we use for Zt in Eq. |^. 

A detailed analysis of K(t) will be given elsewhere. 
Here we note that both parts of Z t {uj) in Eq. ^| have 
the same form. For Z 2 b, the corner frequency uj c — 
l/(Ceff K 2^2D) satisfies Tilu c k B T and we may use 
the kernel of Ref. |l3], while for Z^c we find that 
u> c = VQ/rgCei 2 usually satisfies hui c <C k B T and 
requires different treatment. Since K(t) is linear in 
Re[Z t (cj)], we may calculate K(t) and P(E) sepa- 
rately for and Zg and find the total Ptot{E) as 
a convolution. For Z 2 d, then, we have P2d{E) — 
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where g = Rq/Z 2 y> and B(x, y) is the beta function Hl3| 
For Zrc, we find that 

Knc(t) = - 2 -^Lk B T\t\/h + I [cot(-j^) 
g rc I 2 2k B I t 
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is valid for 2 \ 9 ^ T ~ ""/^' wner e g rc — Rq/{ t iP) and 
t = Rqc^Ij 'V6- From this we calculate 
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complete gamma function. 

To proceed we need an accurate model of Zt(uj). For 
the V g = 7-V curve to be linear at V as 8 /xV, 
Zt/RQ must be nonnegligible at frequencies of order 
eVy/i « 2 GHz. Since then Z 2D « = 20 O, ^ must 
dominate Z t for small ^2D- We therefore consider the 
structure of our leads, which vary in width w from 100 nm 
to 20 /dm. The 100 nm section has length I = 1 //m, 
contributes only 50 17 to Z((0), and is not considered 
further. For the remaining sections with w = 0.4, 1.0, 
10 and 20 /Ltm, and £ = 9, 57, 253 and 375 /im we use 

w R sq /{w + 5.8/i) and q w eeo(w/h + 1.393) where 
/i = 50 nm is the 2DEG depth and e — 13 the dielec- 
tric constant of GaAs to calculate re = 29, 16, 1.9 and 
1.0 Mfi/m and a = 1.1, 2.5, 23 and 46 nF/m. These 
four sections form a cascaded RC line which determines 
Zi{uo). The total Z t (u>) calculated from Eq. || is shown for 
different values of Z 2 d = 1/(4G 2 d) for the pool geometry 
in Fig. [I The cascaded form for Zg(io) is quite complex. 
For our calculations we take Re{Zt(u)] = J2iM Z Rc\ 




FIG. 3: Calculated Z t (u) for (bottom to top) Z 2 n = 258, 
1291 and 4302 Q, for three different models: cascaded RC 
lines (solid), approximate series combination (dashed) and 
Z2U only (dotted). The dash-dotted line is the impedance of 
an infinite line with re = 2.9 Mfi/m and ct = 23 nH/m. 



where the Z^ c are the impedances of the individual sec- 
tions, a very good approximation to the more exact re- 
sult, as shown. Note that this model predicts a signifi- 
cant Z t at 2 GHz dominated by Zi (Z 2 ^) for small (large) 
Z 2 \). For the stripe, Z 2 \^ approaches i?Qpc/2 at zero fre- 
quency and the much lower stripe resistance i? s t r ~ 200 f2 
at frequencies above l/i?QpcCstr where G s t r ~ 0.3 pF is 
its capacitance to ground. At high frequencies, then, Z t 
in the stripe is always dominated by Zi, even for large 
negative V g . We have also shown the impedance for Z 2 d 
alone, and for an infinite RC line with w — 10 fim and 
re chosen to give the correct Ze(0) if the line were finite. 
The latter two models give a small Z t for small Z 2 d at 
the relevant frequencies, and cannot explain the linear 
region in our V g = I-V characteristics. 

Using P 2 d{E) and Prc(E) above, we numerically con- 
volve the P^'(E) for junction j and section i to find 
P { / ] (E) = P[ j) (E) * P 2 U) (E) * P 3 (3) (E) * P[ 3) (E) . We then 

calculate P$(E) = P^ } (E) * P 2D (E) for different G 2D 
and set up a master equation using the rates in Eq. |l| to 
calculate the S-SET current and conductance Gg ET . The 
results for G| ET in the pool geometry are shown as the 
solid line in Fig. ||(a); we scale G§ ET to match the max- 
imum measured value Ggff w 0.34G at G^g x w 6.5G 
but use no other variable parameters. Gg ET agrees rea- 
sonably well with Gset for G2D < G^f x although it rises 
less steeply with G2D- In this regime P 2 d{E) is broad 
and inelastic transitions suppress the coherent supercur- 
rent. For G2D > G 2 I Jf x , G| ET gradually saturates at 
0.47G . For G 2D » G^ x , P 2D (E) « 6(E) [i. e., only 
elastic transitions are likely) and Pi(E) dominates the 
I-V characteristic. No nonmonotonic behavior occurs 
in Gg ET , in agreement with Wilhelm et al. The drop in 
Gset for G2D > G!, 1 ^* must arise from other physics. 
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FIG. 4: Measured (a) and calculated (b) I-V characteristics 
for an unconfined 2DEG (o), V g = -0.3 V (□), and Z 2D = 
1613 (+), 2151 (o) and 6453 fi (x). To fit the data at V g = 
and —0.3 V, we use a c h = 0.07 and 0.05e, respectively. For 
the remaining curves we take a c h = 0. 

We can model the drop by assuming that the measured 
I and Gset are actually averaged over charge states close 
to n g , due to charge motion in the substrate (li). We 
expect charge averaging to be most pronounced when the 
2DEG is unconfined, and least for small G 2 d- Assuming a 
Gaussian distribution of charge states with mean n g and 
variance a c h , we calculate the average (Gg ET ) . We do not 
know the absolute size of a c h, so we assume for the pool 
that a ch = for G 2D < Ggf x . For G 2D > G 2 1 g x we find 
the values of a c h which give (G SET ) = Gset for the pool 
and stripe and plot the results in the inset to Fig. ||(a). 
In both cases a c h is just below 0.04e near V g = V and 
drops near V g — —0.28 V. For the stripe, a c h saturates 
at just above 0.02e, about half the drop for the pool. The 
model seems reasonable given the small a c h required to 
explain the discrepancies with the environmental theory. 

We gain further confidence in the model by comparing 
measured and calculated I-V characteristics for the pool, 
as shown in Fig. [| for SI. For increasing confinement / 
first rises at all voltages, with little or no broadening of 
the linear region (V g = and —0.3 V and Z 2 r> = 1613 fi). 
This corresponds to a reduction in a c h with little change 
in Pt t(E); inelastic transitions in the leads dominate the 
energy exchange. Eventually, a c h — and Z 2 d is large 
enough to affect P to t (E) , causing I to decrease (especially 
at low bias) and the linear region to broaden. The level of 
agreement between the shape and evolution of the mea- 
sured and calculated curves is surprisingly good, given 
the uncertainties involved. While the calculated current 
is much larger than is measured, such discrepancies are 



common in small tunnel junction systems [14 ]. 

In conclusion, we have measured the effects of dissipa- 
tion on transport in an S-SET for which the environ- 
ment can be varied locally. We find good agreement 
with a model in which fluctuations in the leads and low- 
frequency switching between charge states dominate for 
low confinement (large G 2 d), while for strong confine- 
ment (small G 2E ) fluctuations coupled via the capaci- 
tance G 2 d dominate. The model accounts well for the 
evolution of Gset and the I-V curves as G 2E is varied. 
We believe a convolved Ptot = Pi* P2D is likely required 
to interpret the results of the Berkeley group, which may 
explain the discrepancies between their results and the 
scaling theory of Wilhelm, et al. 
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